Effect of Exclusion of Double Occupancies in t-J Model: Extension of Gutzwiller 

Approximation 
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A new type of analytic estimation of the effect of strong 
correlation is developed for the two-dimensional t-J model. 
It is based on the Gutzwiller approximation which gives 
the renormalization of parameters, t and J, due to the 
Gutzwiller's projection operator excluding the double occu- 
pancies. The finite-range correlations are taken into account 
compared with the conventional Gutzwiller approximation 
where only the on-site expectation values are considered. It 
is shown that the essential point of the renormalization is its 
nonlinear dependence on the expectation values of Cooper 
pairs and antiferromagnetic moment. In particular the renor- 
malization factor of J becomes anisotropic in the presence 
of antiferromagnetic moment, i.e., that for the 2-component 
is enhanced compared with that for the xy-component. The 
physical origin of this enhancement is identified as the sur- 
rounding antiferromagnetic correlations of a bond. The self- 
consistency equations for the uniform variational wave func- 
tions are derived and solved numerically. Our result gives 
a reasonable estimate of antiferromagnetic order parameters 
near half filling, in contrast to the conventional slave-boson 
mean-field theory and the original Gutzwiller approximation. 
It is also found that, at half filling, the renormalization of J 
represents some of the quantum fluctuations of the Heisen- 
berg spin system in a different manner from the spin-wave 
theory. For finite doping, our results have some similarity 
to SO(5) symmetric theory. Applications to the inhomogc- 
neous systems such as the vortex state, around nonmagnetic 
impurities, and stripe state are discussed. 

71.10.Fd, 74.20.Mn, 74.25.Dw 



I. INTRODUCTION 

For understanding the basic physics of high-T c super- 
conductivity, the two-dimensional t-J model has been ex- 
tensively studied as one of the most promising and simpli- 
fied models which describe charge and spin dynamics in 
the Cu02 plane. [001 Mean-field theories ||-|(| and nu- 
merical calculations" 7|~|To| have shown that the d x 2_ y 2 
wave superconductivity (SC) takes place in a reasonable 
parameter region in the phase diagram (J/t = 0.3 and 
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the doping rate 6 < 0.3). Therefore, as far as the d-wave 
SC is concerned, the experimentally observed phase dia- 
gram of high-T c superconductors and the results in the 
t-J model agree quite well. 

Recently, however, there are some experiments which 
indicate the importance of the interplay between the an- 
tiferromagnetism and d-wave SC. The most interesting 
phenomenon is a realization of stripe state in some ma- 
terials, which has charge ordering as well as incommensu- 
rate antiferromagnetic (AF) spin ordering. ]TT|— |T^| Other 
potential problems are the AF state induced around non- 
magnetic impurities and a possibility of AF vortex cores. 

In order to study these problems, it is necessary to 
develop a theory in which the d-wave SC and the AF 
correlation are treated in a reliable way. Usual mean- 
field theories generally overestimate the AF long-range 
order, so that they give unphysical results even in the 
uniform case. In this paper we develop a new type of 
analytic theory which gives a reliable estimate of the AF 
correlation as well as the d-wave SC. Our scheme is an 
extension of the Gutzwiller approximation. [5|Jl^| 

In the t-J model the double occupancy of up- and 
down-spin electrons at any site is prohibited. To study 
its ground state, it is natural to use a projected mean- 
field wave function in which the double occupancies are 
excluded. Q| As the mean- field wave functions, a SDW 
mean-field solution at half filling, BCS wave func- 
tion s, MM and a coexistent state of the AF and d-wave 
SC J20|-[22|] have been used. Variational energies are cal- 
culated in the variational Monte Carlo (VMC) simula- 
tions, which treat exactly the constraints of no double 
occupancies. Among the above trial states, the coexis- 
tent state has the lowerest variational energy in the dop- 
ing region 6 < 0.1. @-|| 

On the other hand, the analytic theories in which the 
constraint is treated approximately give large discrepan- 
cies in the estimation of the AF correlation. For example, 
in the slave-boson mean-field theory which takes account 
of the constraint as an average, the AF order is too over- 
estimated and it extends up to unphysical doping rates 
(5 < 0.2). |Q Thus in the slave-boson theory, we are un- 
able to discuss the stripe state, for example, because it 
is stabilized near 6 — 0.125 doping where the slave-boson 
theory gives the AF state even for the bulk. 

The Gutzwiller approximation is more advanced ap- 
proximation than the simple mean-field theories in treat- 
ing the constraint. Renormalizations of expectation val- 
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ues are introduced by comparing the statistical weight- 
ing factors in the wave functions with and without pro- 
jection. ||,^4| As a result, the parameters t and J are 
rcnormalized to gtt and g s J. It has been shown that 
the Gutzwiller approximation gives a fairly reliable esti- 
mation of the variational energy for the pure d-wave SC 
state when it is compared with the VMC results. PJTo] 
However it was shown that there is no region in the phase 
diagram where the AF state is stabilized. [|| This con- 
tradicts with the VMC results. 

In order to clarify the origin of these discrepancies, 
we investigated the VMC data |2^] and found that the 
Gutzwiller approximation has to be modified when the 
AF correlations are present. Based on these observa- 
tions, we extend the Gutzwiller approximation in this 
paper and derive an analytic formalism for the renormal- 
ization factors which reproduces the VMC results. We 
show that it is important to take account of the longer- 
range correlations for the weighting factors, in contrast to 
the previous approximation where only the site-diagonal 
expectation values such as (n ir7 ) are considered. As a 
result, the renormalization factors for t and J have non- 
linear dependences on the expectation values of Cooper 
pairs and AF moment. We think that this is the essence 
of the strong correlation in the sense that the renormal- 
ization appears solely from the projection operators. 

The merit of the present scheme is that it can be easily 
applied to the inhomogeneous systems, where the VMC 
simulations are not so feasible. Thus a reliable analytic 
formulation can be given to the problems, such as the 
stripe state, vortex cores and states around impurities, 
where the interplay between the AF and d-wave SC plays 
an important role. Preliminary results of these problems 
have been described elsewhere. [ ^6|Jl7| , ^5[| 

The outline of this paper is as follows: In Sec. II we 
briefly review the formulation of the original Gutzwiller 
approximation in order to make the present paper self- 
contained. We also show our final results for the extended 
Gutzwiller approximation before going into the details of 
the derivation. In Sec. Ill, the original Gutzwiller ap- 
proximation is extended in order to include longer-range 
correlations. We divide the whole system into cells and 
introduce the configurations in each cell for evaluating 
the weighting factors. Using the general formulation de- 
veloped in Sec. Ill, we obtain the renormalization factors 
for the half-filled case in Sec. IV. We approximate the 
weighting factors in a perturbative way with respect to 
the nearest-neighbor correlations. In IV. D, the physical 
origin of the enhancement of the renormalization factor 
for the ^-component of exchange interaction is discussed 
in a viewpoint of statistical weights of real-space spin 
configurations. In Sec. V, the results at half filling are 
generalized to the less-than-half-filled case. The present 
method is applied to the projected variational states in 
Sec. VI. The self-consistency equations are derived and 
solved numerically to show that they give reasonable es- 
timate of AF long-range order near half filling. Finally 
section VII is devoted to a summary and discussions on 



related problems. 

II. T-J MODEL AND THE ORIGINAL 
GUTZWILLER APPROXIMATION 

We consider the Gutzwiller approximation for the two- 
dimensional t-J model on a square lattice: 

H = -tJ2 Pa{clc ja + h.c.)P G + jJ2 S i- S J> (2- 1 ) 

(ij)cr (ij) 

where {ij} represents the sum over the nearest- neighbor 
sites and Si = c\ a {^cr) a fjCip. The Gutzwiller's projec- 
tion operator Pq is defined as Pq — Ylj(^ ~ 

For this Hamiltonian, we assume the projected BCS- 
SDW mean-field wave function 

|V,)=F G |^o(ALAr f ,M V )), (2.2) 

where A^, A^ f and /x v are the variational parameters 
relating to d-wave SC, AF and chemical potential, re- 
spectively, and ^(A^, A^f , /i V )) is a Hartree-Fock type 
wave function with d-wave SC and AF orders. The wave 
function is a natural generalization of the RVB state 
proposed by Anderson, [jlj It has been shown that the 
coexistent state of AF and d-wave SC has the best varia- 
tional energy in the VMC simulations for the doping rate 

6 < o.i. {jjg-H 

In evaluating the variational energy, the projection op- 
erator makes difficulties for an analytic approach. The 
Gutzwiller approximation was developed for this pur- 
pose. In this method, the effect of the projection is taken 
into account by renormalizations of expectation values as 
follows: 

(4t c jv) = #t(cj CT c itr }o, (Si ■ Sj) = g s {S t ■ Sj) , (2.3) 

where (• • -)o represents the expectation value in terms of 
\ipo) = l^o (A^ AXr, ^ v )), and (• • •) represents the nor- 
malized expectation value in \ip) = Pq\iI>)q\ 

[ ' ~ <V#) (MPgPg\^) ' 1 ' 1 

The coefficients gt and g s are the renormalization factors 
of the expectation values, which we call as Gutzwiller 
factors in the following. Using these notations, the vari- 
ational energy E val = (H) is rewritten as 

£var = (Hcs)o, (2.5) 

where the parameters t and J in H. are replaced with 

t c s = gtt : Jos = g s J- (2.6) 

A systematic estimation of gt and g s was developed 
by Ogawa et al. |is|| whose method is briefly reviewed in 
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II. A. For the case with AF order parameter, it can be 
shown that 



fit 



25(1-5) 
1 - 5 2 + Am 2 ' 



9s 



4(1 - 5) 2 



(1 - 5 2 + 4m 2 ) 2 ' 



(2.7) 



where 5 is the hole concentration, 8 = 1 — n, and m is 
the expectation value of AF order parameter in terms of 
\ipo) defined as 



(-iy 



(<«iT)o - (n 3 i)o) 



(2.8) 



These Gutzwiller factors, however, do not reproduce the 
VMC results, as mentioned in I. We find that this dis- 
crepancy is due to the fact that only the site-diagonal 
expe ctat ion values are taken into account in obtaining 
Eq. (2.7). In this paper we extend the method by Ogawa 
et al. systematically to include the longer range correla- 
tions. We show that it is important to include the effects 
of the nearest-neighbor expectation values, such as 



(2.9) 



in the estimation of gt and g s . This corresponds to taking 
account of the effect of surroundings of the corresponding 
bond. 

Furthermore we find that the Gutzwiller factors g s for 
the exchange interaction have different values for the z 
component (denoted as gf ) and xy component (gx Y )> 



i.e. 



(SiSj)— gf{S*Sj}o, 



(2.10) 



instead of ( |2.3| ) in the presence of m. Actually from the 
VMC calculations, ]2^] we have estimated the behaviors 
of gf and gf Y as a function of m, using the relations 



fJ, 



(SfSJ) 



fJ. 



XY 



( s t s 7)o 



(2.11) 



where the numerators (SfSj) and (S^Sj } are evaluated 
numerically in VMC simulations. It was found that gf is 
enhanced compared with g^ Y and that this enhancement 
is essential for the stabilization of the AF order. In the 
present extended Gutzwiller approximation, this result is 
reproduced. Physically the weighting factor of gf for a 
specific configuration of a bond is enhanced due to the 
effect of the surrounding AF correlations. 

Before going into details of the derivation, we show our 
final results. They are summarized as follows: 



XY = I 2(1 -g) > 7 
1 l-<5 2 + 4m 2 ' 



(2.12) 



and 



„XY 



1 



4m 2 + X 2 



X 2 +4m 2 {l 



6X 2 (l-8) 2 ; 
1 - S 2 + 4m 2 



(2.13) 



where 
a= 1 



4V 



(1 - 5 2 - 

X= 25 2 (A 2 -x 2 
X 2 =2( x 2 + A 2 ). 



Am 2 ) 2 
+ 8mV 



A 2 ) + 4( x 2 + A 2 ) 2 , 

(2.14) 



The above expressions sufficiently reproduce the results 
obtained in the VMC simulations. However we find that 
the slight difference can be improved by using 



X = 2<5 2 (A 2 



X 2 ) + 8m 2 ( X 2 + A 2 ) + 2( X 2 + A 2 ) 2 , 

(2.15) 



instead of Eq. (2.14). A typical m-dependence of gf and 
gf Y are shown in Fig. 1 for the half-filled case, 5 = 0, and 
with A being fixed at the two values, A = 0.02 and 0.18. 
The behaviors in Fig. 1 agree with those obt ained in the 
VMC simulations using the relation in Eq. ( 2.11 ). The 
meanings of the rather complicated Gutzwiller factors 
become apparent in the following sections. 



g s z (A=0.02,m) 



8=0 



g/(A=0.18,m) 




FIG. 



1. Typical m-dependences of the Gutzwiller's renor- 



malization factors, gf and gf Y for the half-filled case, 
(g* Y ) is the Gutzwiller factor for the z- (xy-) component 
of the exchange interaction, respectively. A is the near- 
est-neighbor expectation value of (cj^ct^o in the wave func- 
tion without the projection. It is fixed at values A = 0.02 and 
0.18. m is the expectation value of l/2(n^ — n^) without the 
projection. 

The Gutzwiller factor for the hopping term is given by 



25(1-5) (1 + <5) 2 -4to 2 



2V, 



1 - 5 2 + Am 2 



(1 + 5) 2 - Am 2 



(2.16) 



Figure 2 shows the m-dependences of gf Y , gf and gt for 
8 = 0.12, which can be compared with the VMC results 



3 



in Ref. p2fl. Note here that the slave-boson mean-field 
theory simply gives gt = {b\bj) = 6 and g s = 1. The non- 
linear dependences of , gf and gt on A, \ an d m are 
beyond the slave boson theory. In the following subsec- 
tion we review the original formulation of the Gutzwiller 
approximation which is useful for the generalization in 
the later sections. 



C*2 



5=0.12 



g/(A=0.18, m) 




0.1 



0.2 



FIG. 2. Typical m-dependences of gf , g* Y and gt/5 for the 
doping rate 8 = 0.12 in the t-J model. The meanings of the 
parameters are the same as in Fig. 1. 



A. Original formulation of Gutzwiller approximation 

Let us start with the formulation of the Gutzwiller 
approximation developed by Ogawa et al |l8| when it is 
applied to the t-J model. In this subsection we do not 
consider the AF order. The ultimate aim is to evaluate 
the variational energy 



(H) 



(2.17) 



using \ip) = Pc\4>o)- First, Ogawa et al rewrote the de- 
nominator as (Pq = Pq) 

= i^o\Y[(i - n n ri n )\ij ) 

3 
3 

+ (l-n iT )(l-%;)}|^ ). (2.18) 



The expansion of the product in (2.18) leads to various 
real-space configurations. Therefore {tp\tjj) can be rewrit- 
ten as 



config. j€A j' £Z3 

x II (! -«j"t)( 1 -^"i)l^o), 
3"es 



(2.19) 



where A (B) is a subset of the lattice sites which is singly 
occupied by an up-spin (down-spin) electron and £ is a 
subset which is vacant. The summation is over all the 
possible configurations, i.e., all the possible combinations 
of {A,B,S}. 

In principle, Eq. ( 2.19| ) can be calculated using Wick's 
theorem since \ipo) is a mean-field wave function, al- 
though it is very complicated. The simplest Gutzwiller 
approximation is to estimate (2.19) by using only the 
site-diagonal expectation values, which gives 

config. j£A j'£l3 

x II <(1-%"T)(1- Vl))o- (2-20) 
j"e£ 

with (■ • -)o meaning the expectation value in I'i/'o)- 

For example, the contribution from the site j G A is 
evaluated as 

u A = (n iT (l - n J - i )> = - (1 - -), (2.21) 

where n is the average electron density, n = N e /N, with 
N (N e ) being the total number of sites (total electron 
number). We call lo a as the weight for j site belonging 
to the subset A. Since the hole density is S = 1 — n, we 
obtain 



UJA= UJ B = 



U)E- 



n / n\ 
si 1 "*) 



l-S 2 



2J 4 

Then the expectation value ( 2.20| ) is rewritten as 



(2.22) 



config. 



,N A , Mb. ,Ne _ 



N\ 



, Na. ,N b N 

N A \N B \N E \ L ° A UJb " e 



(2.23) 



where Na is the number of sites in the subset A, and so 
on, (Na + Nb + Ne = N). In this simple case, Na is 
equal to the number of up-spin electrons, so that 



N A = N B = ^, 



Nf 



N -N 



(2.24) 



A similar procedure can be carried out for the estima- 
tion of the numerator in E val . For the exchange term, 
we have 

(iP \P G S+S-p G \^ Q ) 

+(l-n J - T )(l-n ii )}|Vto> 

config. j£-A' j'BB' 

x II (i-Vr)(i-Vi)l*). (2.25) 
fee' 
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Here the summation is over all the possible configurations 
{A! , B',£'} in whic h the two sites £ and to are excluded. 
The evaluation of (2.25) with the site-diagonal expecta- 
tion values for j ^ £,m leads to 



(N — 2)! jf. 



N A ,\N B ,\N E ,\" A ~ l -■'< 



(2.26) 



Since the sites £ and m contain one up-spin electron and 
one down-spin electron, we have Na' = 2 — 1> N B i = 
^ - 1, and N E > = N E = N - N e . 

By combining the estimation of the denominator in Eq. 
( 2.23 ), we obtain the original Gutzwiller approximation 



(S^~S m )- 



(1 + 5) 2 
This gives gf Y = 4/(1 + 5) 2 



(MPaStS m P G \^) 

A 



(2.27) 



B. Antiferromagnetic case 

Ogawa et al JlJ] also considered the above formalism 
in the case with AF long-range order for the Hubbard 
model. The application to the t-J model was carried 
out by Zhang et al. \ 5| In this subsection we follow their 
methods to obtain the Gutzwiller approximation in the 
t-J model. 

In the AF case, the sub-lattices 1 and 2 are distin- 
guished and their magnetizations are defined as 



1 



((%t)o - (^ii)o) = m, 



(2.28) 



for the sublattice 1 and — m for the sublattice 2. Thus 
we denote 



( n jVo= g +m = r > 

Ti 

(%l)o= --m = w, 



(2.29) 



for the sublattice 1, where r (w) means the average elec- 
tron density with the right (wrong) spin direction on the 
sublattice. r and w are exchanged for the sublattice 2. 

Using these notations, the denominator of E vai is eval- 
uated as 



config. 



, N A 2, Nb2, ,N E 2 

UJ A2 W B2 L0 E2 



;f)' 



^ N Al \N B1 \N E1 \ N A2 \N B2 \N E2 \ 

JVai,jvbi 



v , .Nai, ,N B i, ,N E 1, ,N A2 , ,N B 2, ,N E 2 
XU1 A1 UJ m UJ m UJ A2 UJ B2 LO E2 , 



(2.30) 



where the configurations in the sublattice 1 are spec- 
ified by (N A i, N B i, N E i), and so o n. Contrary to 
the previous subsection, Eq. ( 2.30| ) has a summa- 
tion over possible values of Nai and Nm because 
there are only four constraints between six numbers, 
(N A1 ,N m ,N E i, N A2 , N B2 ,N E2 ): 



N A1 



N B1 
N A1 



T' 



N 



A2~ 



2 



N A2 + N B2 
N B i + N B2 



N 

N E2 = -, 

El 
2 ' 



(2.31) 

The weights in the site-diagonal expectation values are 

wai= (%t)o(l - nji) = r(l - w), 
Vbi= (rij i)o(1 - n n ) = w(l - r), 
uj e1 = (1 - n n ) (l - n n ) = (1 - r)(l - w), (2.32) 

for the sublattice 1, and 



WA2 = 



UB1, 
0JA1, 



(2.33) 



for the sublattice 2. 

Since N A i, N B i and TV are large numbers, Eq. 
further approximated by the largest term in the summa- 
tion. By taki ng the partial derivative of the logarithm of 



each term in ( [2.30 ) with respect to Nai, we obtain 



lnAT A1 +ln(^ 



Nai) + IncJAi - lncj A2 = 0, 



(2.34) 



with Nai being the value of Nai wh ich gi ves the largest 
term in the summation. By solving ( 2.34 ), Nai is given 
by 

r(l -w) N e _ r(l - w) N e 



Nai = 



r) 2 



r(l — w) + w(l 
Similarly we obtain 

w(l - r) N e 



2rw 2 



N m 



(2.35) 



(2.36) 



n — 2rw 2 

The summation in ( 2.30) ) is approximated by the term 
with Nai and N B i- 

Again the numerator such as (ip^PcSt S^Pq^q) is 
approximated in the similar way. As a result, the 
Gutzwiller approximation in the presence of AF order, 
to, becomes 



(S^S m ) — 



(MPaSj S~P G \^) 



Nai N B1 _j _ l7 + . 



2 2 
^2 



(n — 2rw) 2 
4(1 - 5) 2 



(1-<5 2 +4to 2 ) 



(Sg~S m )o 

(S 1 / S m )o, 



(2.37) 
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where a useful relation 



2(1 -6) 



n — 2rw 1 — S 2 + 4m 2 



has been used. Equation (2.37) gives 

4(1 ~S) 2 
(1 - 5 2 + 4m 2 ) 2 ' 



XY 



9. 



which is simply Eq. ( |2.7| ). 

It is straightforward to see 



Z XY 

9 s = 9 S i 



and 



9i 



26(1 - 5) 
1 - S 2 + 4m 2 ' 



(2.38) 
(2.39) 

(2.40) 
(2.41) 



III. EXTENSION OF THE GUTZWILLER 
APPROXIMATION: FORMULATION 

The approximation in the preceding section is re- 
stricted to the site-diagonal (on-site) expectation values. 
It was shown that this approximation is not good enough 
to reproduce the AF order obtained in the VMC simu- 
lation. H A straightforward extension of the previous 
Gutzwiller approximation is to consider longer-range ex- 
pectation values. The nearest-neighbor expectation was 
considered before. (2^^^] In this section, we develop a 
general formalism to take account of the longer-range ef- 
fects systematically. 

In evaluating the denominator of E vax , i.e., 
(V'o|-Pg-Pg|' ( /'o} j various configurations have been classi- 
fied depending on the state (|, \ or a hole) on each site 
as {.4, B, £} in the previous section. Instead of this we 
divide the whole system into cells consisting of N c sites. 
(Later we take N c as a large enough number.) For each 
cell, there are 3 Nc configurations because each site has 
1,1 or a hole. We denote these configurations as states 
of a cell. Using these states, the denominator of E vaY can 
be approximated as 



K 



{mp g p g \^)= IK 

all the possible states i—1 

(JL)\ K 

v IK 



{Ni} Ili=l(-^i)' i=l 



(3.1) 



with K = 3^. Here N is the number of cells in the 
i-th state, and u>i is t he we ight of the i-th state cell in 
the analogy with Eq. ( 2.30 ). Explicitly the weight LUi is 
defined as 



3 j' 

x ri(i - ^j"t)(i - "i"i.)i^o>, 



(3.2) 



with j,j' and j" being the sites in the cell. To calculate 
uii is the most important part of the present theory which 
we carry out later. 



The last summation in (3.1) is over all the possible 
values of Ni under the following constraints: 

^n^Ni = —, ^n ii N i = —, 



K 



s Tn hi N i = N-N e 



(3.3) 



where n-j-j («j.i) represents the number of up-spins (down- 
spins) in the i-th state cell, and rthi represents the number 
of holes, respectively. Since a cell has N c sites, a relation 



+ + nu = N c 



(3.4) 



holds. 

Under the abo ve constraints, we look for the largest 
term in Eq. ( |3.l|) as was done in the original Gutzwiller 
appr oximation (see Eq. (2.34)). Since the constraints 
( |3.3|) are rather complicated, it is useful to introduce 
the Lagrange multipliers (/if, A). Taking the partial 
derivative of 



In 



(&)' 



K 

K 

-VlC^n^Ni 

i=l 

we obtain 



iV, 



N e 



K 

i=l 
K 



Ne 

2 



-X^nuN-N + Ne), (3.5) 



N = cj 4 exp(-^|n|i - nin^ - Xrihi)- (3.6) 

/ij = [i, and using the 



From sy mm etry we find /Ltf 
relation (3.4) we get 



Ni = une 



-/iiV e „(/i-A)nhi 



(3.7) 



Furthermore we introduce new variables W and p as 



N 1 

N~r.W' 



Using these variables, N is rewritten as 
N c W 



(3.8) 



(3.9) 



The varia bles W and p are to be determined from the 
constraints (3.2) which are equivalent to 



K _ N K — 

T, N * = w> T, n ^ = N 

8=1 C 1=1 

These conditions become 



N K 



(3.10) 



G 



K 



K 



(3.11) 



It is convenient to classify all the possible states of a cell 
into subgroups which contain j holes (j = 0, 1, 2, • • ■ , N c ). 
(We call these subgroups as j-hole sectors.) Then the 
quantity 



E 



(3.12) 



i with j holes 



represents the total weight of the states in this subgroup. 
Using Wj , we can rewrite the constraints as 



N, 



•sr^ Wj } 



3=0 



J=0 



(3.13) 



The first equality is rewritten as W = J2j Wj%P so that 
W represents the total weight. In the following sections, 
we evaluate Wj and then determine W and p from Eq. 

The numerator in E ml is evaluated in a similar way. 
Let us consider the expectation value of an operator O 
(such as S~[Sm) which operates on a certain part inside 
a cell of the io-th state. (We call this cell as the central 
cell.) Then all the configurations in the other cells are 
classified by {N^} which is the number of cells of the i-th 
state. Then the expectation value becomes 



(iPo\P G OP G \^ ) 



i {N'.} Yii=l(^iY- 



■ N 



K 



IK*<s>< 



(3.14) 



in the analogy with Eq. (2.26). Here (O)i indicates 



the expectation value of O together with terms such as 
Y[i ^i|(l — J^ij.) inside the central cell of the io-th state. 

Determining the largest terms on the right-hand side 
and taking the ratio to the denominator, we finally obtain 



(O) 



(MPgPg\4>o) ' 



E 



w 



(Oh 



(3.15) 



The evaluation of the largest term and the derivation of 
the above result is shown in Appendix A. 

If we consider only the site- diagonal expectation val- 
ues to estimate u>i, Eq. ( 3.15 ) reproduces the original 
Gutzwiller approximation shown in the previous section. 
This is summarized in Appendix B. 



IV. HALF-FILLED CASE 

To obtain the Gutzwiller factors formulated in the last 
section, it is necessary to evaluate the weights Ui for 
each state. Although these weights are generally compli- 
cated, we estimate them assuming that the corrections 
to the original Gutzwiller approximation (reproduced in 
Appendix B) are small. To be more specific, we evalu- 
ate Wj in the lowest orders with respect to the intersite 
correlations \ and A defined in Eq. ( |2.9| ) . Typical values 
for x an d A are less than 0.2, so that the perturbation 
with respect to them will be justified. 

In this section we calculate the Gutzwiller approxima- 
tion for the half filling, since it is less complicated than 
the doped case. All the cells do not contain holes so that 
only the weight, W , is to be calculated. We call this sub- 
set of states as 0-hole sector. After the half-filled case, 
it is rather straightforward to extend the results to the 
doped case. 



A. Evaluation of Wo 

First we calculate the weight u>i of the i-th. state which 
has up-spin electrons, nj,j down-spin electrons: 

w< = (JJ "Jtl 1 II %'i( 1 -% , t))o- (4- 1 ) 

jeA j'eB 

In the lowest order of \ an d A or in the zeroth order, we 
have to take only the site-diagonal expectation values. 
This gives 



wf 3 = [r(l - w)] n '^ HI - r)]"— , 



(4.2) 



where n r ; g ht (n wmn g) means the number of sites on which 
the right (wrong) spices of spin direction is located de- 
pending on sublattice 1 and 2. The summation over all 
the possible configurations gives 

Wo 0) = K 1 - «0 + w (! - r )} Nc = (« - 2rw) N °. (4.3) 

In the next order with respect to \ an d A, we have to 
consider contributions from expectation values of bonds 
in the cell as shown in Fig. 3. For example, consider a 
bond connecting the sites i and j (see Fig. 3(a)) and a 
state in which two up-spin electrons occupy both i and j 
sites. The expectation value in this bond gives 

P TT = ((1 - n a )n lT n iT (l - n iX ))o 

= rw(l — r)(l — w) — rw\ 2 — (1 — r)(l — w)x 2 
-r(l - r)A 2 - w(l - w)A 2 + ( x 2 + A 2 ) 2 , (4.4) 

where we have used the Wick's theorem. Generally \ 
and A can be complex numbers and x 2 , A 2 in the above 
expression mean |x| 2 , |A| 2 implicitly. In the similar way 
we calculate P^i, P^ , P±±. Writing that the left-spin in 
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the subscript indicates the spin on the sublattice 1, they 
are calculated as 



Pi i = Pi 



P 



U 
\ I'- 



ll, 



Pi 



r 2 (l -wf +2r(l - w)x 2 
+r 2 A 2 + (l-w) 2 A 2 
L1 =w 2 {l-r) 2 + 2w(l-r) X 2 
+w 2 A 2 + (1 - r) 2 A 2 



(x 2 
(x 2 



A 2 ) 2 , 
A 2 ) 2 . 



(4.5) 



J □ 



(a) 



(b) 



(c) 



FIG. 3. Real-space diagrams for the evaluation of the 
weighting factors of the 0-hole cells, (a) Contribution from 
a bond which gives X in Wq 1 ' and (b)(c) neglected diagrams. 

In the estimation of I^i 1 ' we need their summation 



P TT + P n + P iT 



P 



1 1 



(n - 2rw) 2 + X, 



(4.6) 



with X defined in Eq. (2.14), where we have used the 
relations r = n/2 + m and w = n/2 — m. The first term 
(n — 2rw) 2 on the r.h.s. is the zeroth order contribution 

which has been included in Wq '. In this sense, X is a 
kind of connected (or cummulant) contribution coming 
from the real-space diagram in Fig. 3(a). The actual val- 
ues of X is roughly 1/20 so that our perturbation scheme 
can be justified. Denoting the number of bonds in a cell 
as Nb, we have a contribution to Wq as 



W (1) = N b X(n - 2rw) N "- 2 , 



(4.7) 



where the factor (n — 2rw) Nc ~ 2 comes from the contri- 
butions of the other N c — 2 sites in the cell except for i 
and j. 

In the similar way, we evaluate the higher order con- 
tributions to Wo. By considering two bonds in the cell 
we approximate their contribution as 



Let us discuss the neglected contributions to Wq. Fig- 
ures 3(b) and 3(c) show real-space diagrams whose con- 
nected expectation values can contribute to Wo- Al- 
though their contributions are not included, i t ca n be 
shown that they are smaller than the terms in (4.9). For 
example, the contribution from Fig. 3(b) is in the order 
of Nb6 2 x A and iVf,m 2 x 4 etc., so that it is smaller than 
Wq 1 ^. The contribution from Fig. 3(c) is in the order of 
X 2 Nb which is smaller than fl4.8| ) by a factor 1/Nb- 

There is ano ther effect neglected in Wo'- The higher 
order terms in (|4.9| ) become less and less correct because 
the number of the bonds giving X is not exactly N b Cj 
due to their overlapping. However the error is again in 



the order smaller than N h Cj. In this sense, Eq. (4.9) is 



a kind of a summation of most dominant terms in the 
order of Xm 3 b . 



B. Evaluation of gf Y 

In order to obtain the Gutzwiller factor for S^S^ at 
half filling we need to calculate 



(4.10) 



according to the general formulation ( 3.1 5| ) . Note that 
W = Wq and the central cell (which includes sites I and 
m) does not contain holes, i.e., rihi = 0. The average in 
a central cell of the i -th state, {S^S~)i , is evaluated 
in a similar perturbation scheme as in Wo- In the lowest 
order of \ an d A, we have (Fig. 4(a)) 



{SjS m ) {n - 2rw 



,N c -2 



(4.11) 



where the factor (n — 2rw) Na ~ 2 comes from the contri- 
butions from the sites in the central cell other than the 
sites I and m. 



m 



m 



W (2) = Nb C 2 X 2 (n-2rw) 



(4.8) 



where N h C% = Nf,\/(Nb — 2)!2! is the number of choices for 
the positions of the two bonds in a cell. The summation 
of these series leads to 



J2 N b C j Xj{n-2rw) N °-V 

3=0 



(n — 2rw) 



X 



(n — 2rwY 



(4.9) 



(a) 



(b) 



FIG. 4. Real-space diagrams for the evaluation of (S~[ S^)- 
(a) In the lowest order calculation, we have (Sf S^o for the 
sites I and m. (b) In the next order, the contributions of X 
come from the bonds which are not connected to the sites I 
and m. The dotted lines and the solid line are the excluded 
bonds. 

In the next order of x an d A, there are contributions 
of X from the bonds in the cell. This leads to 
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N b X(S+S-) (n - 2rw) 



(4.12) 



where N b is the number of bonds which are not connected 
to the sites I and m directly (Fig. 4(b)). In general N b 
depends on the shape of the cell. However, if we choose 
the large enough cell, we have 



N b = N b - 7, 



(4.13) 



as is evident from Fig. 4(b). 

Although the dotted bonds in Fig. 4(b) do not give X, 
they may still give different contributions to (S^~ 5*"}^, 
such as 



(1 



*m' J. 



(4.14) 



where ml is a nearest neighbor site of m and (■ ■ -) c means 
a connected expectation value excluding the disconnected 
terms such as (iS , /S'~)o(n m '|(l — n m 'i))o- (The discon- 
nected expectation values have been already taken into 
account in Eq. (4.11).) By calculating (4.14), however, 
we find that there are no connected conributions: Ac- 
tually the two electron operators in Sg~ = cLc^ have to 
make contractions with the electron operators on the site 
m in the Wick's expansion, which gives just the discon- 
nected expectation value. 



By summing up the series of corrections such as ( 4.12 ) 
we obtain 



'y* J (Sg' s m )i 

= (S+S m ) {(n-2rw) N c- 2 

+N b X(n - 2rw) Nc - 4 + ■ ■ ■} 

= (S+S m ) (n-2rw) N c- 2 (l + 



(n — 2ru>) : 



N,, 



(4.15) 



Combining with Wo, the final result for the extended 
Gutzwiller approximation for S 1 " becomes 



(Sc S m ) 



(Sg~S m )o 
(n — 2rw) 2 



1 



X 



(n — 2™) s 



-(N b -N b ) 



(4.16) 



Since the Gutzwiller factor is defined as the ratio between 
(SfS~) and (S/S~)o, we have 



1 



(n — 2rw) 2 
a -(N b -N b ) 

(n — 2rw) 2 



X 



(n — 2rw)' 



~(N b -N b ) 



where we have put 



a = 1 



X 



(n — 2rw) 2 



(4.17) 



(4.18) 



which appears frequently in the following discussion. If 
we choose a large enough cell like A c >> 2, we obtain 



XY 

9s = 



(n — 2rw) 2 



(4.19) 



according to Eq. ( 4.13 ). Physically the factor a~ 7 with 
a > 1 represents the exclusion effect of the bonds as 
shown in Fig. 4(b). 



C. Evaluation of gf 

The Gutzwiller factor for the z-component of Sg ■ S m 
is calculated in the similar way. The lowest orders with 
respe ct to \ an d A give the similar results to ( 4.11 ) 
and (4.12), in which (S^S~)o is replaced with (S , |S , ^ l ) . 



However, for (S^S^)^, an additional contribution ap- 
pears from the diagram in Fig. 5(a) which did not give 
contributions in {S^S~)i . This is the reason why the 
anisotropy for the Gutzwiller factor between xy compo- 
nent and z component appears, which is the most impor- 
tant feature of our extended Gutzwiller approximation. 

Assuming that the site I and ml in Fig. 5(a) are in the 
sublattice 1, we have a contribution 

(SfSninm'tQ - n m >i) +n m 'i[l - n m < T )}) c 

= {Sj;)o(S* n n m ^(l - ri m >i) + S^n m 'i(l - n m / T )) c 



:(P- 



rr 



P 



n 



Pi 



It 



Ai) 



= -m 2 X 2 , 



(4.20) 



with Xi defined in Eq. ( 2.14 ). Since the number of bonds 
connected to the site m is 



No 



N b — N b — l 



(4.21) 



the contribution to JA {^e^m)io turns out to be 

-2N 2 m 2 X 2 (n-2rw) N < ! - 3 : (4.22) 

where the contributions from the bonds connected to the 
site £ are also included. 



Ill 



m 



I 



m 



m 



(a) 



(b) 



FIG. 5. Real-space diagrams contributing to {SjtSm). A 
kind of connected expectation values of (a) and (b) give ad- 
ditional contributions to gf, compared with g* Y . 
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Similarly the diagram in Fig. 5(b) gives a contribution 

(5fn*/ t (l -n t >i) + Sfn ti {\ - n^ T )) c 



-m 2 J^. 



(4.23) 



By counting the number of possible combination of the 
bonds, we have 



N£m 2 X£(n - 2rw) 



N c -4 



(4.24) 



The higher order terms of are calculated as before, 
which become 



1 



(5f S* > (n - 2rw) N °- 2 a N » 

-2N 2 m 2 X 2 (n - 2rw) N "~ 3 a N '' 
-N%m 2 Xl(n - 2rw) N '-*a<), (4.25) 



with TVj and NfJ being the numbers of the bonds which 
are not connected directly to the diagrams in Figs. 5(a) 
and 5(b), respectively. They are shown in Figs. 6(a) and 
(b). In the large enough cell, we have 



JV£= N b - 10, 



Ng= N b - 13. 
Combining with Wq, we obtain 



a *>) 



(4.26) 



(n — 2rw) 2 



2N2 m 2 X 2 a-^-^) 



n — 2rw 



(n — Lrwy 



(4.27) 



m 



m 



m 



(a) 



/ 

—4 

V 



m 



(b) 



FIG. 6. The bonds excluded for the contribution of X, 
which are connected directly to £, m and m! sites in Fig. 5(a) 
and £,m,m' and £' sites in Fig. 5(b), respectively. 



In this case, (S^S^) is not proportional to (SfS^o. 
Since the Gutzwiller factor is the ratio between the two, 
we have some nontrivial contribution from the second 



and third terms on the r.h.s of ( 4.27 ). Using 



(4.28) 



the final expression for the Gutzwiller factor turns out to 
be 



Z XY 

9s = 9s 



1 



Am 2 + X 2 



X 2 +4m 2 {l + 



N 2 X 2 



n — 2rw 
where we have assumed 

N b - Ng = 2(N b - Nl). 



(4.29) 



(4.30) 



For a large enough cell we have -/V2 = 3 and Nb — N^ = 3, 
so that 



Z „XY 

9s =9 S 



1 



Am 2 



X 2 



X 2 +4m 2 {l 



3X> 



n — 2rw 



(4.31) 



The above formula is one of the most important results 
in this paper. Typical m-dependences have been shown 
in Fig. 1. Let us here check some limiting cases. When 
m = 0, we have 



= gf = 4a- 7 , 



(4.32) 



which reproduces the original Gutzwiller approximation. 
For small m, we obtain 



gf = <?f y h + 4m 2 



6 



9X 2 



n — 2rw 



(n — 2rw) 2 



(4.33) 



It is apparent that gf has an enhancement as a function 
of m, which gives the reasonable AF long-range order at 
half filling. 



D. Physical meaning of the enhancement of g s 

In the usual interpretation of the Gutzwiller approx- 
imation, a comparison is made between the probabili- 
ties of spin configurations in the wave functions with and 
without the projection operators. When we con- 

sider only two sites I and m for Si ■ S mi the summation 
of the probabilities of spin configurations | TT)> I TDJ IT) 
and I 11} in the wave function without the projection is 
given by P^-f + + Pj t + Pg which was calculated as 
(n — 2rw) 2 + X in Eq. (4.6). On the other hand, in the 



wave function with the projection, we have always one of 
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the above four states at half filling, and thus the proba- 
bility is equal to 1. The ratio of these probabilities gives 
the Gutzwiller factor 



-XY 



flf = 



1 



(n - 2rw) 2 +X' 



(4.34) 



However this Gutzwiller approximation does not give a 
reasonable answer compared with the VMC results as 
mentioned before. 

For the physical understanding of the results in the 
previous subsections, it is necessary to go beyond the 
two sites, £ and m, because the enhancement of gf comes 
from the diagrams in Fig. 5. For this purpose we consider 
the spin configurations on three sites as shown in Fig. 7. 
To calculate the probabilities of these configurations is 
straightforward. For example for Fig. 7(c) we have 



= 2r(l - w)P n - 2r 3 (l - w) 3 



(4.35) 



hances the probability of this configuration, 



where the second term on the r.h.s. appears to avoid the 
double counting of the zeroth order with respect to x an d 
A. We can see that the presence of the third site ml en- 

f 

.T I, 

because P^ > Pjj = Pjj > P^ in the presence of AF 
order. As a result, the weight of the configuration f J, 
on t and m sites is increased in the wave function. This 
causes the enhancement of gf . 



m 



The denominator becomes 

(n - 2rw){(n - 2rw) 2 + 2X}, 
while the numerator becomes 



(4.37) 



Pi 



11/ 



(n-2rw)(P n -P n -P iT 
+2m(P TT - P n + P iT - P u ) + 4m 2 (n - 2rw) 
= 4(n - 2rw){S z t S z m ) - 4m 2 X 2 . (4.38) 

Althou gh th e denominator is larger than the two-site case 
in Eq. (4.34), the second term in the numerator enhances 
{SfS* n ), which is mainly from the contribution of the 
configuration 7(c). Using 4(5|5^) = -(4m 2 + X 2 ) and 
taking the ratio, we obtain 



(SfS^) (n - 
X 2 



1 



1 



2rw) 2 + 
f 4m 2 (l 



2X Am? 
X 2 



2rw 



X 2 



(4.39) 



This is essentially the result obtained in the previous sub- 
section for gf . Since we have co nsid ered only the site m! 
as the third term, N 2 in Eq. ( 4.2E ) is replaced with 1 
in this case. The first factor l/{(n — 2rw) 2 + 2X} cor- 
responds the exclusion effect which was represented by 
a -{N b -N b ) m gXY _ p rom t ms analysis it is apparent that 



the enhancement of gf as a function of m is due to the 
increase of the probability of || in the presence of the 
AF circumstances. 



Z A— A A— t A-^ 



(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



(g) 



(h) 

FIG. 7. All the spin configurations on £, m and m' sites, 
which contribute to gf . The probability of the configuration 
(c) is the largest because of the AF correlations (£ and m! 
sites being on the sublattice 1). The effect of the surrounding 
AF correlations (in this case on the site m!) is the physical 
origin of the enhancement of gf . 

To understand this effect more quantitatively, we cal- 
culate (S^S^) directly from the configurations in Fig. 7 
as follows: 



(a) + (b)-(c)-(d)-(e)-(f) + (g) + (h) 
(a) + (b) + (c) + (d) + (e) + (f) + (g) + (h)- 



(4.36) 



V. LESS-THAN-HALF-FILLED CASE 

A. Evaluation of the total weight W 

In the presence of holes, the Gutzwiller factors are cal- 
culated similarly as in the half-filled case. In this subsec- 
tion we obtain the total weight, W. Since the weight for 
the zero-hole sector Wq has been calculated, we have to 
evaluate 



Wi 



E 



UJ, 



(5.1) 



■ith 1 hole 



and W 2 and so on. 

In the lowest order of \ an d A, we have 



W[ 0> = N c (l - r)(l - w)(n - 2rw) 



N, 



(5.2) 



where the factor N c comes from the possible position of 
the hole in a cell and the factor (1 — r)(l — w) comes from 
the expectation value of the hole position, ((1 — n^)(l — 
nj|)) . In the next order, by counting the number of 
bonds which contribute X, we obtain 



N c (l - r)(l - w)N lb X{n - 2rw) 



(5.3) 
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where N± b means the number of bonds which are not 
connected the hole site as shown in Fig. 8(a). In the 
large enough cell we have 



N lb = N b - 4. 



(5.4) 



(a) 



(b) 



FIG. 8. The bonds excluded for the contribution of X, 
which are connected directly to (a) the hole site, and (b) the 
sites of two holes. 

Although the dotted bonds in Fig. 8(a) do not give 
X, they give alternative contributions to W\ in the same 
order. This contribution is calculated as 

((1 - n iT )(l - n u ){% T (l - nji) + n n {l - % T )}} 
= (n-2rw)(l-r)(l-w)+Y, (5.5) 

with 

Y = 6(1 + S)( x 2 - A 2 ) - 4m 2 ( x 2 + A 2 ) - 2( x 2 + A 2 ) 2 . 

(5.6) 

Here i represents the site of the hole and j is one of the 
nearest neighbor sites. Since the number of the dotted 
bonds in Fig. 8(a) is (Nb — Nib), the additional contribu- 
tion to Wi becomes 



N c (N b - N lb )Y{n - 2rw) 



AT -2 



The higher order terms with respect to X can be taken 
into account in the similar way as before and we obtain 



W t = N c (l - r)(l - to) (n - 2rw) 



+N c {N b - N lb )Y(n - 2rw) N ^ 2 a Nh 



(5.8) 



where N b appears because the diagram in Fig. 8(a) ex- 
cludes the same number of bonds as in Fig. 5(a). We 
assume that 



N b = N lb - 3, 

which leads to 

Wi= N C U1 - r)(l - w) + 
x(n~2rw) N ^ 1 a Nu . 



(5.9) 



*»-*»y fl -3 
n — 2rw 



(5.10) 



Strictly speaking, when the hole is located on the bound- 
ary of the cell, the weight should be different. However 
we make an approximation of large enough cell. 



Continuing the similar arguments, we obtain 
W 2 = Nc C 2 |(l - r)(l -w) + Nb ~ Nlb Ya- 3 



n — 2rw 



x(n-2rw) N "- 2 a N2 



(5.11) 



for the two-hole sector, where N21, is the number of bonds 
which are not connected to the two sites of holes (Fig. 
8(b)). We have neglected the case in which t wo ho les 
come to the nearest neighbor sites. From Eq. ( |5.11 ) it 
is apparent that the j-hole sector, Wj, can be generally 
approximated as 



Wj = Nc C 3 z 3 {n - 2rw) N °- j a N i-» , (5.12) 



with 



z=(l-r)(l-w) + Nb f lb Ya- 3 . (5.13) 
n — 2rw 

Finally approximating that Nj }b = N b — 4 j, the con- 
straints which determine W and p (Eq. J3.13] )) become 

E ^P 3 = i(n - 2rw +pza-yca N > = 1, 



3=0 



3=0 



w 



w 



pza (n — 2rw + pza ) 



4\N -l„N b 



= SN C . 

From these constraints we obtain 

6(n — 2rw) 4 
P= a , 



(5.14) 



W= 



n — 2rw 



(5.15) 



(5.7) for less-than-half-hlling. 



B. Evaluation of g* Y and gf 

In the similar way we estimate gf Y . In the case when 
there arc j holes in the central cell, we have 

io with j holes 

= Nc -2 C jZ i(n - 2rw) N °- 2 -ia">-V(S+S-) . (5.16) 

By summing up all the sectors with different nu mber of 
holes, we have (according to the general formula ( 3.15 )), 



JV -2 



{Sg~S m )— E^ — (SfS m )i 

j—0 io with j holes 



1 

w 



(n - 2rw + pza-^^a^iS+S^o 



W- N >HStS-) . (5.17) 



n — 2rw 
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From this we obtain the Gutzwiller factor 

2 



I) 



n 



2rw 



a 



-7 



(5.18) 



This is a simple g ener alization of the results at half filling 
obtained in Eq. ( hlG ) . 

Although the estimation of gf is essentially the same as 
in the half-filled case, there are some additional diagrams 
to be taken into account which are shown in Fig. 9. The 
detailed calculations are summarized in Appendix C. The 
final result is 



z XY 
9 S = 9 S 



1 



Am 2 + X 2 



r NoXoTl , 

4m 2 {l + (1 

n — 2rw 



3^ 2 



r\ -31 



(5.19) 



The quantity p is given in Eq. (|5.15| ). For gf, we use the 
lowest order approximation 



5(n — 2rw) 



2S, 



(5.20) 



since n — 2rw = 1/2 and z — 1/4 in the lowest order with 
respect to 8, \ and A. This gives our final expression for 



at, 



given in Eq. ( 2.13; ). 
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Qm' 

lb 

(f) 



FIG. 9. Real-space diagrams which give additional contri- 
butions to gf in the presence of holes. The detailed calcu- 
lation for each configuration is summarized in Appendix C. 



where the expectation value of the correlated hopping, 
T, is calculated as 

T= ((1 - nii)c\ r Cjf(l -%0)o 
= (l-r)(l- W ) X -X 3 -xA 2 



(5.23) 

It is easy to show that for the j + 1-hole sector 

N ^ 2 CjTa N ^ b (n - 2rw) N "- 2 - j z j . (5.24) 
Therefore summing up all the sectors, we obtain 

JV -2 



3 W>=W E N c -2C j Ta N ^(n-2rw) N °- 2 -3zipj+ 1 



3=0 

n 



n — 2rw 



Ta ~(N b -N b ) 



p. 



(5.25) 



Since (cj«Cj-[-)o = x, the Gutzwiller factor g t — 
( c !t c jt)/( c 1t c Jt)o becomes 



9t= 



{(i- r )(i- w )-^}o-(^-^)p 
n — 2nt; y 2 

n <5{(l-r)(l- W )-f } (Nh _j} h)+4 



n — 2rw 



(5.26) 



where we have substituted the value of p obtained in Eq. 
( PI ). The diagrams in Fig. 5 give contributions of the 
order of S 2 y 2 or 8m 2 x 2 so that they are neglected. 

Finally comparing the definitions of X and Y, we use 
an approximation 



Y 



X 



(5.27) 



for the small S case. Using this approximation, the quan- 
tity z can be rewritten as 





z= (1 


-r)(l 




C. Evaluation of gt 


= (1 


-r)(l 


-w){l- 


As for the Gutzwiller factor for the hopping term, g t , 


= (1 


-r)(l 




we have to calculate 


= (1 


-r)(l 


— wja 



N b - N- 



lb 



(1 — r)(l — w)(n — 2rw) 



Ya 



(5.28) 



((1 -n a )ct T c jT (l -n jU ))i 



(5.21) 



Since the hopping term needs at least one hole, the lowest 
sector is the one-hole sector. For this sector, we have 

2q with 1 hole 

= ((1 - n 4i )4 T c jT (l - % i )) a^(n - 2rw) N ^ 2 

= Ta^ b (n-2rw) N "- 2 , (5.22) 



Here we are considering 5, m and x as small quantities. 
Substituting this approximate z into gt we have 



<5n (1 — r)(l — iu) 

.9< = 7 w r^ -0 

n — 2rw (1 — r)(l — w) 



X2 

_2_ n -(N b -N b )+8 



25{1-S) (1 + 5) 2 - 4m 2 - 2A 2 



l-5 2 +Am 2 (l + <5) 2 -4m 2 



a. (5.29) 



This is our final expression for g t given in Eq. (2.1q). 
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VI. OPTIMIZED VARIATIONAL STATE 



and the exchange term is replaced with 



In this section we calculate the energy of the varia- 
tional state using the Gutzwiller factors obtained in the 
previous sections. We will show that the variational en- 
ergies and the magnitudes of order parameters agree with 
the results in VMC simulations. 

In the variational state l^} = Pg\^o(A^ , A^ f , /x v )) in 
Eq. (Q, |V>o(Ay, A^ f ,/i v )) is a Hartree-Fock type wave 
function with the d-wave SC and AF orders. It is ex- 
pressed as Q 

imam,m v )h n (4 s) +4 5) 4? t ^il)io) 

fe,«(=±) 



= n4 ex p 



k , s 



V V k ^)t.,(«)t 

_.(«) a fcT a ~ki 



k,s u k 



where 

„,(±) 



dVk 



,(±) 



El 



{±E k - (j> 



Ay f 2 , 



|0), (6.1) 



(6.2) 
(6.3) 



(±E k - + (A y %)2 



(-k = — fryfc, Ik = 2{cosk x + cosky) and r\ k = 2(cosfc,; - 
cos ky). The annihilation operators d k s ^ are related to the 
electron operators through the following unitary transfor- 
mation, 



(6.4) 




)■ 


f Ot k< j 


-Pk* \ 


f CAka \ 




I (3k* 


a k cr J 


I CBka J 



with 



Q-ka 
Pka 



(6.5) 



Here CAkaic-Bkcr) are annihilation operators of an electron 
on the A(B)-sublattice and a represent T(+l) and !(-!)• 
The wave vector k is limited to half of the Brillouin zone 
where e& < 0. We can confirm that \ipo) is a vacuum of 
the annihilation operators which diagonalize 

E [E Uk(c\ ka c B ka + h.c.) 

k a 

-(/i V + CJ^ t )c\ k(y CAka ~ (M V " <rA%f)CBka C Bk<r} 



"AX?7fe(c j 4_ fei C BfeT + C B ~klCAk] + h.c.) 



(6.6) 



JSi-Sj = gf Y J(S*S 



SVS*) 



g?JS?S* 



(6.8) 



where gt , gf Y and gf are the Gutzwiller factors obtained 
in the previous sections. When we apply the mean-field 
theo ry t o 7i c ff, we obtain the similar Hamiltonian as in 
Eq. (|6.6|) but with the replacements 



t 



t e ff 



Aj^ J off A v , 



with 



2J c z s m Y , 



T z - n z T 
• J eff — 9s J - 



(6.9) 



(6.10) 



(6.11) 



Note that in the usual mean-field theory the self- 



consistency equations give x 



X 



A and 



m — m, where Xi A and m are the expectation 
values A = (c^cj^o, \ = ( c L c j<?)o, and m = 
j(— 1)* ((n?f)o — (fhl)o) used in the previous sections. 
We will see shortly that x V and X e tc are slightly differ- 
ent due to the dependence of the Gutzwiller factors on 
X, A and m. 

Using the wave function |-0o) , the expectation values 
become 



-E 

8N tX 

-E 

SN tX 

-E 

2A ^ 



F k r\k 



v /(±£ fc -M v ) 2 +^ 



(*eff + J cS X V H(±E k - M V ) 



, ±E k ^{±E k -^f+F } 



2J z B m v (±E k - M v ) 



k,± 



±E ky /(±E k -n 



V\2 



F 2 

r k 



1 E 



(±E k 



y/(±E k -^)2 +F * 



with 



F k = J cS A rjk, 



E k = V(*eff + JeffX V ) 2 7fc 



(2Jf ff m 



V\2 



(6.12) 
(6.13) 
(6.14) 
(6.15) 



(6.16) 
(6.17) 



where N is the number of sites and the summation over 
k is limited to half of the Brillouin zone. 
Using these expectation values we obtain 



In order to clarify the correspondence to the mean-field 
theory, let us consider the effective Hamiltonian, 7i e s in 



Eq. (2.5). In H c s the parameter t in TL is replaced with 
t C B = gtt, (6.7) 



E vm - — (Heff)o 

= -8Nt cSX -4,N J ofr (A 2 + x 2 
-2NJ z s m 2 - Li{Nn - N c ). 



(6.18) 



14 



By taking the derivatives with respect to the variational 
parameters, A v ,x V , and to v , we find the following self- 
consistency equations: 



1 



dH c ff 
\ a a / 



8NJ eS x dA 

1 



m 



v 



TO 



/I 



8iVJ eff 
1 



<9x 



>o, 



<9to ' 



1 



TV on 



(6.19) 



where the partial derivative of H c fi is applied to the 
Gutzwiller factors, g* Y ,gf and g t . 

Figure 10 shows the self-consistent parameters A, x 
and m satisfying ( |6.19| ) as a function of the doping 
5 = 1 — n for J/t = 0.3. Note that these are the ex- 
pectation values in the wave function without the 
projection. The dashed line in Fig. 10 represents the re- 
sults when the AF order is suppressed, i.e., to is fixed to 
zero. It is apparent that the presence of the AF order 
has a small effect on the expectation values A and \. 
The coexistent state between the d-wave SC and AF or- 
der parameters is stabilize up to the doping rate S = 0.1. 
This is consistent with the VMC simulations. 




FIG. 10. The self-consistent parameters A, x and m as a 
function of the doping rate S = 1— n for J/t = 0.3. The dashed 
line represents the results when the AF order is suppressed 
i.e, m is fixed to zero. 

The actual expectation values in the wave function 
PgI^o) with projection are different from A,x and to. 
For example, the expectation value of cj^Cja- is 



(6.20) 



In the similar way, the actual expectation values in 
Pa\ipo) are obtained using corresponding Gutzwiller fac- 
tors. For m cxp we repeat the similar arguments as for gf 
to obtain 



TO exp = - (nit _ n ii) 



2rw 



g m m, 

{N b -N lb )X 2 n { 
n — 2rw 



(6.21) 



For the expectation value of cLcL, we have to take care of 
the fact that the number of holes changes in the cell. Here 
we use the average number of holes as an approximation. 
In the same level of approximations as for g t we obtain 

A ft xn = (c. 



j i 
dn 



9A-- 



(1 



r) 2 + (l- W ) 2 +X 2 (Nh _X h)+8 



n — 2rw 2(1 — r)(l — w) 

26(1-5) (1 + S) 2 + Am 2 + 2X 2 
1-6 2 + Am 2 (1 + (5) 2 - 4to 2 



(6.22) 



Figure 11 shows the actual expectation values, A e , 
and TO eX p as a function of the doping 5 = 1 — n for J/t 
0.3. 



0.4 



m exv =g m m 




A e xp=,?AA 



A A A A 



••I*!! 



0.2 



and m e 

in the optimized (or self-consistent) wave function with the 
projection operator as a function of the doping rate 5 = 1 — n 
for J/t = 0.3. 



VII. SUMMARY AND DISCUSSION 

In this paper we have developed a new approach 
for studying the effect of strong correlation or the 
Gutzwiller's projection in the two-dimensional t-J model. 
It is based on the extended Gutzwiller approximation, in 
which the effects of longer-range correlations are taken 
into account. These correlations play important roles for 
the interplay between the AF and d-wave SC. Let us sum- 
marize our main results and discuss related problems. 

(1) Generally the expectation values with respect to 
the projected wave function are strongly renormalized 



15 



due to the exclusion of double occupancies. In the slave- 
boson mean-field theory, this effect is taken into account 
by assuming the replacement t — > 6t and J — > J. How- 
ever we have shown that the renormalization is not so 
simple. First of all the renormalization factor for the ex- 
change interaction is anisotropic, i.e., = g* Y J ^ 
Jjj = gf J in the presence of AF moment. Further- 
more gt , gf Y and gf have nonlinear dependences on the 
expectation values A, \ and to. In this sense the ex- 
tended Gutzwiller approximation is beyond the simple 
slave-boson mean-field theory. We think that essence of 
the strong correlation is contained in these Gutzwiller 
factors, because they stem solely from the projection. 

The physical meanings of the Gutzwiller factors are 
clarified: 

(1) The factor a~ 7 appearing in g^ Y and gf repre- 
sents the exclusion effect shown in Fig. 4(b). When we 
calculate (SiSj), the surrounding six bonds cannot make 
contribution of nearest neighbor correlation X. As is ev- 
ident in Fig. 1, the factor a~ 7 reduces the value of J e g 
even for the case m = 0. 

(ii) The enhancement of gf as a function of m is the 
most important feature. This was observed in the VMC 
simulations. p2| Here we have identified its origin as the 
effect of surrounding AF correlations by discussing the 
probabilities of spin configurations in Fig. 7. The en- 
hancement of gf is due to the increase of fj, configuration 
caused by the AF circumstances. 

(2) We have studied the projected variational states 
as an application of the new Gutzwiller approximation. 
It is shown that our scheme reproduces the results in 
the VMC simulations. ]2^] The enhancement of gf as a 
function of to is essential for reproducing the coexistent 
state of the AF and d-wave SC orders found in VMC. 

mm 

(3) At half filling (<5 = 0), there are several interesting 
features. In the original Gutzwiller approximation, the 
AF state at half filling was not obtained || which was 
unphysical. However in the present Gutzwiller approx- 
imation, the AF state with d-wave SC order parameter 
becomes the most optimized state as obtained in VMC. 
Note that the expectation value of the SC order parame- 
ter, A exp , is zero because of the Gutzwiller factor <?a = 
at half filling, although A and the variational parame- 
ter A^ are nonzero. The self-consistent values which we 
obtain are (Fig. 10) 



This is a reasonable magnitude and close to the Monte 



Carlo result for the Heisenberg model 



0.31 ± 



A v -~ y 
A 



v_ 



X' = 0.080, 
X = 0.163, 
0.017, 
0.147. 



(7.1) 



The actual expectation value with the projection, m Gxp , 
is (see Fig. 11) 



xp = g m m = 0.373. 



(7.2) 



0.02). |2J 

The relations A v — \ V and A = \ are the manifesta- 
tion of the SU(2) symmetry at half filling. j2|,| Actually 
the self-consistency equations for x an d A become the 
same at S — 0. It is worth while noting here that, due to 
the SU(2) symmetry, the coexistent state of the AF and 
d-wave SC is equivalent to the 7r-flux state with AF long 
range order discussed by Hsu. |2(J 

The self-consistent value, m — 0.147 in Eq. (7.1), is 
smaller than the result m — 1/2 in the simple mean- 
field theory. Instead it is slightly larger than the value 
at which the Gutzwiller factor gf has a maximum, as 
shown in Fig. 1. If to is increased from the self-consistent 
value, gf decreases. This indicates that m is determined 
so as to minimize the exchange energy —2Njf s m 2 = 
—2NgfJm 2 by optimizing the energy gain due to the 
long-range order and the energy loss due to the reduc- 
tion of gf. This situation is, in some sense, similar to 
the quantum fluctuations of the Heisenberg spin system 
discussed in the spin-wave theory: The spin fluctuation, 
which reduces m from the mean-field value 1/2, leads 
to the enhancement of the Gutzwiller's renormalization 
factors to gain the energy. 

(4) For less-than-half-filled case, the difference between 
the variational parameter A v and the expectation value 
A exp is remarkable: A^ is finite and moreover it increases 
as the doping rate S decreases, while A exp is proportional 
to 8 near half filling. This is because the Gutzwiller fac- 
tor, <7a, is proportional to S (see section VI). We interpret 
that A v is the BCS-type energy gap observed in scan- 
ning tunnel spectroscopy pp| |34| or in break junctions 
|p5| , because it is the parameter embedded in the wave 
function even at half filling. The excitation spectra will 
have a large energy gap corresponding to A v , while the 
true long-range order, A cxp , is reduced due to the pro- 
jection or the strong correlation. The increase of A^ in 
decreasing S is consistent with the dependence observed 
exp er iment ally. 

(5) We discuss here the relation to the SO (5) theory. 
|E3 In our formulation, the combination to 2 + A 2 appears 
frequently. If the SO (5) symmetry is exact in the t-J 
model, the free energy will have a systematic dependence 
such as 

F(A 2 + to 2 ) or F(A 2 xp + m c 2 xp ). (7.3) 

For example, the numerator in a has a combination 

16(to 2 + A 2 + x 2 )- (7.4) 

Also the factor giving the enhancement of gf contains 

4m 2 + X 2 =4to 2 + 2A 2 + 2x 2 . (7.5) 

Although these combinations remind us of the SO (5) 
symmetry, our Gutzwiller approximation does not show 
exact symmetry. 
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Nevertheless a tendency similar to the SO (5) predic- 
tion can be seen from the Gutzwiller factors in Fig. 1. 
Comparing g* Y and gf for the cases with A = 0.02 and 
0.18, we find that g^ Y and gf are larger for smaller A. 
This is mainly from the exclusion effect of a, because 
even at to = this effect is observed. Therefore if we 
consider a situation where the d-wave SC order parame- 
ter is suppressed, then the Gutzwiller factor is enhanced 
causing the increase of AF moment. This gives the sim- 
ilar phenomena predicted in the SO (5) theory. |l^] 

(6) One advantage of the present theory is that it is 
easily applied to the inhomogeneous cases, such as the 
stripe state, vortex cores, and magnetic states around 
nonmagnetic or magnetic impurities, where the interplay 
between the AF and d-wave SC correlations plays an im- 
portant role. Since our Gutzwiller approximation gives 
a reasonable estimate of the variational energies in the 
presence of AF correlations, a reliable analytic formula- 
tion can be given. 

The simplest way of applying the present scheme to 
inhomogeneous problems is to assume that gf Y ,gf and 
g t for each bond are determined locally from the 

expectation values Aij,Xij an d TOi,TOj for the bond. In 
this case, if the d-wave order parameter is reduced around 
vortex cores, impurities or stripes, then the Gutzwiller 
factors g^ Y and gf are enhanced locally as expected 
from Fig. 1. This effect causes the local development 
of AF correlations, which can be observed experimen- 
tally. From these viewpoint, preliminary calculations for 
the vortex cores [^6[0 and stripe states |2f| have been 
published elsewhere. 

Of course the VMC simulations give more accurate 
evaluation of the variational energies, if they are used 
for the inhomogeneous systems like stripe states. How- 
ever there are some difficulties in applying the VMC. 
Firstly we have to treat a fairly large unit cell to study 
the slowly varying order parameters especially near half 
filling. For example, the incommensurability in the 
stripe state is close to (it, ir) so that the period of the 
stripe pattern becomes fairly long. In this case the VMC 
simulations become difficult. Furthermore the choice of 
the functional form of the trial state in VMC is restricted, 
since only the small number of variational parameters can 
be used practically. On the contrary, our scheme based 
on the mean-field-type Gutzwiller approximation can be 
used in a fairly large system sizes. Moreover the order 
parameters on all the bonds, Aij,x-ij an d m i can be op- 
timized in the similar sense to unrestricted Hartree-Fock 
theory. Therefore we can search for the microscopically 
optimized variational states in our scheme. |^6| , |l7| , p5| 
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APPENDIX A: DERIVATION OF THE 
EXPECTATION VALUE OF AN OPERATOR 



The evaluation of Eq. (3.14) is slightly complicated 



The constraints for {V 4 '} are 

i=l i=l 



i=l 

J2 n hiN-=N -N e -n Uo , 



K 



(Al) 



instead of < J3. 3| ) because N- represents the number of cells 
of the i-th state except for the central cell. 

In the same way as in the denominator, the largest 
term is given by 



— N u. 

W c w 



* AT W KP > ' 



(A2) 



with slightly different values of W and p' because of the 
difference of the constraints. From the constraints (Al), 
we can see that N[ satisfies the relations 

K — N K — 



Therefore if we define the difference ANi = N! — Ni, we 
have important relations 



K 



K 



5^AJVi = -l, J2n hi ANi 



(A4) 



i=l 



By use of these, the r atio b etween the numerator and 
the denominator in Eq. ( 3.15 ) is calculated as follows: 



0-- 



(iJjo\P g OP g \^ ) 
(iPo\PgPg\iPo) 

i \N C ) i=l iv i' i=l 



l 1=1 



{O)io 



K 



V — TTi — — 

^ v Ml n c w 



-ANi 



(O)io 



E-f--V EAJVi xp-s- A ^a) i0 

^ N \N C W I y K ' 



^ w w 



(A5) 



This gives Eq. Q3.15| ) 
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APPENDIX B: REPRODUCTION OF THE 
ORIGINAL GUTZWILLER APPROXIMATION 

In this appendix we show that the original Gutzwiller 
approximation described in II is reproduced if we assume 
only the site-diagonal expectation values in the general- 
ized formulation in section III. 

First we calculate the weight u>i of the i-th. state which 
has 7i|i up-spin electrons, down-spin electrons and 
nhi holes. If we take only the site-diagonal expectation 
values, we have 

Ui = [r(l - w)] n ^ ht [w(l - »•)]*""■« [(1 -r)(l - w)] n ", 

(Bl) 

where n r i g ht (^wrong) means the number of sites where the 
right (wrong) spices of spin direction is located depend- 
ing on the sublattice 1 and 2. Then the total weight in 
the subgroup with j holes (Eq. ( gig )) can be calculated 
exactly as 



E < 

i with j holes 



Ne Cj[T{\ -w) + w(l - r)] N °- j [(1 - r)(l 
Nc Cj(n - 2rw) N °- j (l - r) j (l - w) j , 



(B2) 



where N c Cj is the number of choices of the positions of 

j holes. 

Using these Wj, the constraints ( 3. 13| ) have simple 
forms as 



3 

E^ 



3=0 

N c 



w 



w 



{n - 2rw + (1 - r)(l - w)p} Nc = 1, 



3=0 

x{n- 2rw + (1 - r)(l - w)p} Nc 
These can be solved easily to give 
S(n — 2rw) 



SN C . (B3) 



P 



n(l — r)(l — w) 



W = 



n — 2rw 



(B4) 



Finally the expectation value for St S m , for example, be- 
comes 

(ik\P G StS-P a \1>o) 



{MPgPg\^o) 



N c -2 



w 
w 



3=0 



{Sg S m )i g 

C j (n-2rw) N °- 2 - j (l-r] 



= —{n - 2rw + (1 - r)(l - w)p} N ^ 2 (S+S m ) 



(n — 2rw) 



(B5) 



which is exactly the same as the results in the original 
Gutzwiller approximation in section II. B, Eq. (2.37). 



APPENDIX C: EVALUATION OF Gf FOR THE 
LESS-THAN-HALF-FILLED CASE 

In the similar way to g* Y , we evaluate gf . In the 
zero-hole sector, we have 



E {S z t S* m ) io ={S z t S* rn ) {n-2rw, 

io with holes 



N c -2 a N b 



-2N 2 m 2 X 2 (n - 2rw) N "~ 3 a N "-~ :i 

-N 2 m 2 X 2 (n - 2rw) N ^ i a^- 6 . 

(CI) 

For the one-hole sector, there are five contributions from 
the diagrams in Fig. 9. The first three diagrams are sim- 
ple extensions of the zero-hole sector. However Figs. 9(d) 
and 9(e) are new- type contributions due to the presence 
of a hole. For the diagrams in Figs. 9(d) and 9(e) we 
calculate 



(£^(l-rw T )(l-fw.0) £ 



2 21 



x(l - n m 'i){l - n m /|)) £ 

-m 2 X 2 . 
o 2 



(C2) 



Using these expectation values, we obtain 
E] (^t S m )io 

zo with 1 holes 



-- Na -2 Ci(S z t S z J z{n - 2rw) 1 ^- i a 



N c -3„N lb 



- Na - 3 C 1 2N 2 m 2 X 2 z(n - 2rw) 
-N c -±CiN 2 m 2 X 2 z(n - 2rw) 
+N 2 m 2 X 2 (n - 2rw) 



N c -5 a N lb -6 



N c -3 a N b -3 



*(1 



+N 2 m 2 X 2 (n - 2rw) N ^ 4 a Nb - 6 . (C3) 

In the two hole sector, we need to calculate a diagram 
in Fig. 9(f) which is — im 2 Xf. Taking account of these 
/ diagrams and counting the higher order terms with rc- 

t ™'§p ec t to X, we approximate as 
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iq with 2 holes 

=n c -2 C 2 (S^S z J z 2 (n - 2™)^"^ 
- Nc - 3 C 2 2N 2 m 2 X 2 z 2 (n - 2nu) Ar °- 5 a^- 3 
-N c -iC 2 N%m 2 Xlz 2 {n ~ 2rw) N ^ 6 a N2 "- (i 
+N c - 3 C 1 N 2 m 2 X 2 z(n - 2rw) N ^ i a N ^ 
+ Nc ^ 4 C 1 N 2 m 2 X 2 z{n - 2rw) N ^ 5 a Nlb - 6 



(C4) 



Finally generalization to higher order contributions 
leads to 



1 Nc 



3=0 io with j holes 



N c -2 



= W^ p3 N a -2C j {SfS z m ) z j (n-2rw) N ^a N ^ 



N c -3 



E P j N c -3Cj2N 2 m 2 X 2 z 3 (n-2rw) 

3=0 
V c -4 

]V P 3 N c -iC ] N 2 m 2 X 2 zHn-2rw) N ^ 



II" 
1 



Nc-3-j a N j:b -3 



3=0 
N c -3 



+ W ^ pl+1 ^-sC j N 2 m 2 X 2 z j {n-2 
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+ W ^ p3+1 ^-iO j N 2 m 2 Xiz j (n-2 
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— 2rw 
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(S!S z m } a-( N »- N ^ 
i 2N 2 m 2 X 2 (l-^)a-^-^ 
l N 2 m 2 X 2 (l-^) 2 a-^-^ 



-(N b -N b ) 



X- 



2 2 
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{1 + 
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4 1 n - 2rw y 2 

Thus the Gutzwiller factor is 
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(C5) 



4m 2 + X 2 
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1 n-2rw y 2> ' 
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